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Abstract 

We consider a physical system constituted by a finite chain of point masses con- 
secutively linked by linear springs and dashpots. At one of the end points acts an 
external control force aligned with the chain and the system is observable by the 
position of the other end point. We show that, whatever is the number of the point 
masses, if the sequence of the elastic constants is proportional to the sequence of 
the dashpot constants, then the mechanical system is completely controllable, com- 
pletely observable, completely reachable and completely re-construictible, in the 
sense of control theory. 

Key words: Control, Dashpots, Springs, Mass chain, Point masses chain, 
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1 Introduction 

1.1 Premise 

In [1] a finite chain of mass points consecutively linked by linear springs is 
considered; the mechanical system is externally controlled by a control force 
aligned with the chain, applied at one of the end points, and is observable 
by the position of the other end point. There it is shown that such a system 
is completely controllable, completely observable, completely reachable and 
completely re-construictible, in the sense of control theory, for all choices of 
the physical parameters. 



Here we consider a spring-mass system similar to the one described above but 
includes viscous damping which is caused by the presence of dashpots connect- 
ing consecutive points of the chain. We show that such a system is completely 
controllable, completely observable, completely reachable and completely re- 
construictible provided the sequence of the dashpot constants is proportional 
to the sequence of the elastic constants. This is true for any choice of the 
number of points, of their masses, and of the elastic constants of the springs. 

1.2 On linear systems with control inputs 

Following the same notations of [2], we consider the continuous linear first 
order system 



where z, u and y are real vectors of respective dimensions n, m and p, z = 
dz/dt, F is the (n x n) state matrix, g is the (n x m) input matrix, h is the 
(p x n) output matrix, and d is the (p x m) direct transmission matrix. In 
literature the two equations above are usually referred as state equation and 
output equation, respectively. 

A mechanical system E represented by equations (l)-(2) is said to be ([3]-[5]) 

(A) completely reachable if for each state z there are an instant t > and an 
input function «(.), defined in [0, t], such that the solution z(.) to equation 
(1) satisfies z(0) = 0, z(t) = z; 

(B) completely controllable if for each state z there are an instant t > and 
an input function «(.), defined in [0, t], such that the solution z(.) to equation 
(1) satisfies z(0) = z, z(t) = 0; 

(C) completely observable if every initial state z(0) can be determined by 
observing the future outputs y(s) for s e [0, r] for some r > 0; 

(D) completely re-construictible if every actual state z(t) can be determined 
by observing the past outputs y(s) for s E [t — r, t] for some r > 0. 

By the linearity, for any continuous system (l)-(2), 

(E) the controllability space and reachable space coincide and do not depend 
on the lenght t of the time interval; (F) the observability condition and re- 
construictibility condition are equivalent and do not depend on the lenght r of 
the time interval. 



z = Fz + gu , 
y = hz + du , 



(1) 
(2) 
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Here we consider strict casual systems with one scalar input and one scalar 
output; hence 

d = 0, m =\=p, (3) 



g is a column vector of dimension n, h is a row vector of dimension n, and 
equations (1), (2) become 

z = Fz + gw, (4) 

y = hz. (5) 



We use the following theorem, which holds for system (4)- (5 ), with one scalar 
control input and one scalar output (see [2]). 

Theorem 1.1 The mechanical system S, represented by (4)-(5), is completely 
reachable and completely observable if and only if the polynomials 

det(zI-F) and hadj(zI-F)g (6) 



have no common root. 



1.3 Mechanics of spring-mass chains with viscous damping 



We consider the chain of N > 2 mass points Pi, P 2 ,...,Pn of respective 
masses mi, 7712, . . . , m^. Each point interacts only with its nearest neighbours 
since any two consecutive points Pj, Pj + i are linked by a linear spring hav- 
ing elastic constant kt > and natural lenght £j. Moreover, there is viscous 
damping since consecutive mass points are connected by dashpots of con- 
stants Ci > (see Figure 1 and Figure 2). We denote Zi the displacement of 
the point Pj from its equilibrium position, so that Zi — for each % — 1 to N 
corresponds to no deformations in the springs. Of course, this system has iV 
degrees of freedom and the Zi are free coordinates for it. 

The equations of motion of the mass points P, have the form 
d 2 

mi 175 Zl = kl ( z<2 ~ z i) + c i(zz ~ zi) + ft Xt > ( 7 ) 



d 2 

m i~^£2 Z i = ki—i{Zi—i — Zi) + ki(Zi+i — Zi) 

+Cj_i(ii_i - Zi) + Ci{z i+l - zi) (i = 2, . . . , TV - 1) , 



(8) 
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d 2 

mN ~dfi ZN = ^ n ~ 1 ^ Zn ~ 1 ~ Zn ^ ~^ c n-i{zn-i — z n) + f J 
Here we study 



ext 
N 



(9) 



(i) the controllability of the state system by the input control u = fl xt , the 
external force applied to the first mass point P±, and 

(ii) the observability of each state z by using as output the position y = zn 
of the last mass point. 

Hence we put fff 1 = and by introducing the variables zn+i = -j^Zi, equations 
(7)-(9) give rise to the first order system of 2N linear equations 



d_ 

di 



z i — z N+i , i — 1, . . . , N, 



d 



mx—ZN+i = h(z 2 - Zi) + Ci(zjv + 2 - z N+1 ) + u, 



(10) 



mi ~dt ZN+i ~ ^ i ~ 1 ^ Zi ~ 1 ~ z ^ ^ Zi + 1 ~ Zi ) + 
Ci-i{z N+i -i - z N+i ) + d(z N+i+1 - z N+i ) (i = 2, . . . , N - 1) (12) 



m N —z 2N = k N -i(z N -i - z N ) + c N -i(z 2 n-i - z 2 n) ■ 
at 



(13) 



In order to write equations (10) to (13) in matrix form, we put 



X(xi, x 2 ,---, %n-i) = 



XI 


XI 







m\ 


m\ 




XI 


Xi+X 2 


X2 





m 2 




m 2 





X2 


X2+X3 


X3 


m 3 


rn 3 


m 3 

















V 



o 
o 



XN-2 

XN-2 X N _2+%N-1 



mjv-i 




mjv-i 








m N 





_ £jv_]_ 

m N -i 

XN-l 
m N 



Kjy — X(fci, k 2 , ■ ■ ■ , fcjv-i) ) — X(c 1; c 2 , • • • , cjv-i) , 



(14) 



and 
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'N 



In 

-Cn 



(15) 



where Oat is the square zero matrix of order N and Ijv is the identity 
matrix of order N ; moreover we put 



Z\ 



Z2 



\ Z 2 N ) 



I o ^ 



g 





1/mi 




V 







... 1 ... 



(16) 



where 1/mi is in (N + l)th place of column g, of dimension 2N, and 1 is in 
iVth place of row h, of dimension 2N. 

Hence equations (10)-(13) in matrix form rewrite as: 



z = Fz + git, 
y = h • z . 



(17) 
(18) 



2 Two useful lemmas 



For any square matrix A we use the notation 
= DetA . 



(19) 



In order to prove the two key theorems of the present paper, i.e. Theorem 3.1 
and Theorem 4.1, we need the two lemmas below. 
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Lemma 2.1 For any choice of the positive integers a, (3 and of the real 
number p, let 

O ax p be the a x (5 zero matrix, 

I a be the identity matrix of order a, 

Bp be any matrix of order (5, and 



A 



O a x(3 Pla 
Bp Op xa 



(20) 



Then 



= p°(_l)°(0+2) \\ B p\ 



(21) 



Proof. By expanding A in the cofactors of the 1st row yields 



\\A\\ = p(-l)W 



0( a -l)x/3 Pla-l 
Bp 0/3x(a-l) 



(22) 



Repeating this step to the matrix 

0(q_i) x /3 pla-i 
Bp Op x ( a -i) 



(23) 



consecutively for each % = 2, . . . , a, yields 



\\a\\ = pH-iyxv+y 



0( a ~i)xf3 pla-i 
Bp 0/3x(a-i) 



(24) 



and thus for % = a we obtain equality (21). 



Q.E.D. 
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Lemma 2.2 Let 7 be a positive integer and let A be a square matrix of order 
27 — 1 having the form 



A 



zl 



7-1 



ROW-y^ l(0) 

UpTr(ai, a 2 , ■ ■ ■ , a 7 -i) UpTr(b x , b 2 , ■ ■ ■ , 6 7 -i) CoZ 7 _i(#) 



(25) 



where (i) itow 7 _i(0) is £/ie row of zeros of dimension 7—1, (22) UpTr(xi,x 2 , . . . , x 7 _i) 
denotes any upper triangular matrix having the numbers X\) x 27 • • • , <£ 7 — 1 ^ 
i/ie main diagonal and (in) CoZ 7 _i(#) denotes any column of dimension 
7 — 1 . TTien 



(-l) 7 (ai + 2&i)(a 2 + ^62) • • • (a 7 _i + 2fc 7 _i) 



(26) 



Proof. For each % — 1, . . . , TV — 1 multiply the (N — 1 + i)t/i column of A 
by z and add the result to the ith column; we obtain a matrix B having the 
form 



B = 



-U 



7X(7-1) J 7 

UpTr(aj + zbj : 1 < j < 7 - 1) UpTr(bj : 1 < j < 7 - 1) CoZ 7 _i(#) 



.(27) 



Now for i = 1,...,7, multiplying the ith row of -B by suitable quantities 
and adding the results to the (N + 1 + i)th row, one obtains a matrix B' 
having the form 



B' 



O 



7 x( 7 -l) 



UpTr(dj + zbj : 1 < j < 7 — 1) 0( 7 -i) 



x 7 



(28) 



By the properties of determinants we have \\A\\ = \\B\\ = \\B'\\ . By Lemma 
2.1 with p = — 1 we have 



||B'|| = (_i)7(_!)7(7+i) \\U P Tr( aj + zbj : 1 <j < 7- 1)|| ; 

thus, since 7(7 + 1) is even, 

= ||B'|| = (-1) 7 \\U P Tr( aj + zbj : 1 < j < 7 - 1)|| . 



Now equality (26) holds since the determinant of a upper triangular matrix 
equals the product of its diagonal elements. 
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Q.E.D. 



3 Form of the adjoint polynomial 



Theorem 3.1 In connection with any positive integer N > 2 consider the 
square matrix F N of order 2 N defined by (14)-(15); moreover 

(i) let a N = zI 2 n — F N , thus 



a N 



z I N -I N 
K N C N + zI N 



(29) 



(ii) let 



h N = 



... 1 ... 



(30) 



be the row of dimension 2N , with components 5 Nt , i = 1,2, ...,2N (5 = 
Kronecker's delta); 

(in) let 



9n 



... 



... 

mi 



(31) 



be the column of dimension 2N, with components ^S^^ 1 , i — 1,2, ... , 2N. 



Then 



h N adj(a N ) g N 



(fa + ZCi) (k 2 + zc 2 ) ■ ■ ■ (fcjv-l + zc N ^i) 



mi m 2 . . . rriN 



(32) 



Proof. The polynomial hNadj(aN)gN equals c N N+1 jm\ , where c N N+1 is 
the (TV, N + 1) component of the matrix c = adj(aN); hence we have 



mi h N adj(a N )g N = c N N+1 = (-1) N+ ( N+ V \\A N+1 N \\ = - \\A N+1 N \\ ,(33) 



where A N+ i N is the submatrix of obtained by deleting the (iV + l)throw 
and the Nth column. Now we observe that A N+ i N has the form (25), that 
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IS, 



l N+l N — 



in 



z In-i 
Rown-i(0) 

UpTr(a u a 2 , . . . , a N ^) UpTr(b u b 2 , . . . , 6jv-i) Co^v-i(#) 



,(34) 



with 



ai = - 



61 = - 



m 2 

Cl 

m 2 



a 2 = — 



6 2 = - 



m 3 
c 2 

m 3 ' 



, ajv-i 



, &JV-1 



CjV-l 



(35) 
(36) 



and Co/jv_ 1 (#) is any column of dimension N — 1. Now Lemma 2.2 yields 
HA/v+i TV || = (-1)^(01 + z&i)(a 2 + zb 2 ) ■ ■ ■ iflN—i + zb N -i) , (37) 



and thus equalities (33) to (37) yield 



h N adj(a N )g N 



^ C NN+1 ^ ^ ^ + ^)( a2 j_ zt)^ , , , (Qj\f_i + zfejy-i) 



mi 



mi 



- ( ^AT+l+AT-l + ZC l) ( k 2 + ZC 2 ) ■ ■ ■ (fcjy-1 + gCjy-l) 



mi m 2 . . . m^v 



Q.E.D. 



4 Form of the characteristic polynomial 

Let us remind that the characteristic polynomial Pn(z) of Fn is defined by 
Pn{z) = \\a N \\ , with a N = zI 2N — F N . (38) 

The next theorem characterizes the characteristic polynomials of the mechan- 
ical systems considered here. 

Note that the chain of mass points which is considered here has no springs 
connecting the first point and the latter to a fixed external wall. 

For this reason the characteristic polynomial has the twice root zero. 
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Theorem 4.1 In connection with any positive integer N consider the square 
matrix F N of order 2N defined by (14)-(15), and put 



b N = zcn + kN ■ 



(39) 



Then the equalities 
Pi(z) = z 2 , 



A) = i, 



D N (z) = P N {z) + -JLD N ^{z) , N = 1, 2, ... , 

m N 



(40) 
(41) 



P^W = {z 2 + —)Pn-i(z) + z* ^ 



m N 



m N -i 



D N _ 2 (z), N — 2,3... , (42) 



inductively define the sequence {Pn(z)} n>1 of the characteristic polynomials 
of F N . 



Moreover, let 

{1} N = z 2 I N + zC N + K N + W N , 



(43) 



with 



Wi 



N 



O(AT-l) CoZ(iv_i)(0) 



(44) 



T/ien the quantities D k in (40)- (42) are the determinants of matrices {l} k : 
Dk = \\{l}k\\ ■ (45) 



Proof. We prove the theorem by induction. Hence we firstly note that the 
inductive formulae (40)-(42) hold for iV = 2; in fact, see (15), for 



2 h 
-K 2 -C 2 



(46) 



by (38) we have 



a 2 



zl 2 -I 2 
K 2 (zl 2 + C 2 ) 



(47) 
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and 



P 2 (z) = \\a 2 \ 



(z 2 + z^ + ^ 

v mi mi 



) ( ~ c l _ i-I. 
/ v mi mi 



ci _ fci 



TO2 



) (z 2 + z^ + ^-) 

ni2 ' \ m,2 m2 ' 



mi m 2 I \mi m 2 , 



(48) 



It is easy to verify that using (40)- (42) with N = 2 one arrives at the same 
result. That is, the expression of P 2 (z) given by its definition (38) coincides 
with the expression given by equalities (40)-(42). 

Now, in order to use the induction principle, assume the inductive hypothesis: 
let (40)- (42) hold for an arbitrarly given integer iV; then consider 



Fn+i — 



On+i In+i 
—Kn+i —Cn+i 



(49) 



with Kn+i and Cn+i defined accordingly with (14), (14). 
By (38) we have 



Pn+i(z) — ||djv+i| 



zIn+i —In+i 
Kn+i (zIn+i + Cat+i) 



(50) 



Now, for each i = 1, . . . , N, we multiply the (N + 1 + i)th column of the 
matrix cln+i by z and add the result to the it/i-column; we obtain the matrix 



On+i In+i 
(z 2 In+i + zCn+i + Kn+i) (zIn+i + Cat+i) 



(51) 



Hence for each % = 1, . . . , N + 1 multiplying the ith row of the latter matrix 
by suitable quantities and adding the results to the (N + 1 + i)th row, we 
obtain the matrix 



On+i —In+i 
(z 2 In+i + zCn+i + K N +\) On+i 



(52) 
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Of course, by the determinant properties, 



l N+l 



|OjV+l || — Pn+i{z) 



(53) 



By (52)-(53) Lemma 2.1, for a = f3 = N + 1 and A = ajv+i, yields 



*AT+1 



(_l)JV+l(_l)(^+l)(^+3) 



Pjv+i^) .(54) 



Since JV + 1 + (N + 1)(N + 3) = (N + 1)(N + 4) is even for each N, we have 

(_l)"+l(_l)(*+l)("+3) = 1 5 and thus 



z 2 In+i + zCn+i + K N+ i 



Pn+i(z) 



(55) 



Now, by expanding the matrix z 2 In+i + zCn+i + -Kjv+i in the cofactors of its 
last row, we obtain 



Pn+i(z) = \\z 2 I N+ i + zC N+1 + K N+ i\ 



:-i) 2JV+2 f z 2 + z^^ + J*- ) \\{i} N \\ + (-i) 2N+i 

V mAr + i mjv+i / 



thus, 



+ + -^)I|{1WI + + -r^) ll{2}.|| , (57) 



m N+l m N+1 



m N+1 m N+ i ' 



where, for 



a-,- = z + z— h — — 



mo- 



rn A 



(58) 



the matrices {1}at, {2}at are defined by 



{1} 



N 



ZiZ 2 



{2} 



N 



Z1Z3 



(59) 
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with 



Zi, Z2 and Z% respectively given by 



(z 2 + z^ + ^ 

\ mi mi - 



(_2_ci fci \ 

^ m2 m2 ' 








(_ z _£l «i ■ 

mi mi - 





v mo mo ' 



y m 3 m 3 / 







7712 
_£3 JS3 



7714 



rri4 - 












CN-3 



«AT-3 



mjV-3 "7JV-3- 



y CJV-2 

' mjv_i 







fcjV- 2 

mjv- 



7) 








-2; 



fcjV- 



CiV-2 _ 

mjv_ 2 mjv_ 2 
«AT-2 



I _ fc JV-i ~\ 
mjv mjv ' 







) 

/ Cjy-i fcjy-i \ 

V mjv-i m N -i' 

OlN-1 








fcjV- 



CJV-3 _ 
mjv-3 mjv-3 

ajv-3 

cjv-2 fc 



mjv_i 



mjv- 









CJV-2 



«JV-2 



7) 



,-2 



mjv-2 mjv„ 2 ' 

«AT-2 
Cjv-1 fcjv-i \ 











mjv 



mjv 



( „ cn _ fejv 'j 

\ mjv mjv ' . 



For instance, for Af = 4 we have 



^ 2 + z^ + -^' 

v mi mi • 

v ni2 m2 ' 








i-z^ - ^) 

\ mi mi ' 

( z % _|_ ~ ci+c 2 1 fci+fc 2 ' 
V m2 m2 • 

(-2 -SL - -^) 

v 777,3 m 3 / 
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[z z + z 



m 2 
mz 



+ 



V 014 rru I 



k2 \ 

k2+k-s \ 
mz i 

ki 

rrn i 








~ C3 _ fea ' 

m-3 m,3 • 



> 2 + + 



A:3+fc4 % 

m,4 / 



v ni2 m,2 ' 



/ 2 + ca+ca + 

V m.3 



k2±Ma ) 

mz ' 



' — Z -£2 JS3.) 

\ m4 m,4 / 









(z 2 + z^- + -M - 

mi mi ' v m,i m,i 



ni2 mo / 



mi 
fci_ 



H12 • 








mi mi ' 
Oil 

(-Z&- - 

v m,3 mi / 









' — z .£2 k '2 ) 

mo mo > 



m.2 

a 2 



kz 








v r«4 m4 • 



'_ Z C3_ _ 

mz 



«3 



[z 2 + Z^ + ^' 

v mi mi • 



' ni2 






fcjj 

m2 • 



mi / 



' mi 
, C2 &2 



m3 



m3 



_;r_£2 k 2 \ 

m,2 ni2 ' 






a 2 

v mi ma ' v mj 



(z 2 + 2^ + ^) 
\ mi mi / 








CjV-2 

mjv-2 



m JV _2 - 



"iV-2 
Cjv-1 ' 



mjv 



m/v 











, C N 
' mjv 



fcjy ■ 
mjv - 



14 



[z 2 + z^ + -M 

>> mi mi > 



(69) 



&N-2 \—Z ) 



/ cjv-i 



v mjv 



+ 



fciv-i 



by the linearity property of determinants with respect to columns, we can 
express ||{1}jv|| within (57) as 



II{1} 



Z 2 Ijq + zCjq + K]\[ 



+ \\E N \\ = P N {z) + \\E_ 



N 



(70) 



Again by the linearity property of determinants we have \\En\\ + ||{2} iV || = 
since the last column of the matrix En + {2}^ vanishes; hence, taking into 
account equalities (56) through (68), equality (57) yields 



Pn+i(z) 



.2 ''v k N \ D . 



zT + z- 



+ 

m N +i m N+1 ) 



P N (z) + z 2 \\E 



N 



(71) 



Replacing determinant ||.Ejv|| by its expansion in the cofactors of the nth 
column in equalities (71) and (70), respectively yields 

P N+1 (z)=(z 2 + z-^ + -^_)p N (z) + z 2 (z^L + ^L) {1}^ (72) 



m N+ i m N+ i I 



m N m N 



and 



2 Cat t kN 



:iu=^)+^(v+^) {i>iv-i 



(73) 



where the matrix {1} 7V _ 1 is defined by (43) for % — N — 1. 

Thus we have proved (40) through (42). 

Lastly note that by (59)- (61) equalities (43), (44) yield (45). 



Q.E.D. 
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Remark 4.1 By the above inductive definition, the polynomial Pn(z) can be 
written in the form 



P N (z) = (z 2 + 



ZC N - X + k N -i\ „ , N , 2 ZC N-l + k N -i 



m N 



P N ^(z) + z 2 



(Pn-2(z) + ... 



• • • + (P 2 (*) + (z + • • • )• 

m 2 ™- 



m 3 



mi 



(74) 



Indeed, replacing back with equalities (40)-(41) in (42) yields (74). 

Moreover, the above inductive definition implies the following expression for 
Pn{z) : 



( 



2N-2 

Pn{z) = E E 

p=0 2M+fl=0 
0<M;R<p 



\ 



E 



l<r 1 <r 2 <...<r M <JV 
l<s 1 <s 2 <...<s J j<JV 
\l<tl<t2<...<t M +R<N 



k Tl k r2 . . . k rM c Sl c S2 . . . c SR 
m h m t2 ...m tM+R 



z 2N - p .(75) 



For example, for N = 2 we have P 2 (z) = (z 2 + S£ ^) Pi(z) + z 2 ^^D , 



thus 



P 2 (z)=z* + z 3 (^ + ^)+z 2 + 
\mi m 2 J \mi 



TO 2 , 



For iV = 3 we have 



P S (z) 



z z + 



zc 2 + k 2 ' 
m 3 



zc 2 + A: 2 
m 2 



Pi(z) + 



zci + ki 
mi 



D 



thus 



Ps(z) 



^ 2 _^ zc 2 + k 2 



z* + z*(^ + ^)+z 2 K + h. 



mi m 2 



mi m 2 



2 zc 2 + k 2 I 2 zci + h 
+z z H 



m 2 



mi 



and 



PJ Z ) = Z e +Z 5 C cl + SL. + + ^) +z 4 ( fci. + fcL + fca. + Ja. 

J \ / \mi "12 "*2 "*3 / \™1 "*2 "12 "*3 



miTO2 



m, 2 m 3 , 
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1^3 / fcic 2 I fcic 2 _|_ fc 2 ci _|_ fcic 2 _|_ fc 2 ci _|_ fc 2 ci \ , 2 / fcifc 2 I fcifc 2 I kik 2 
\miBi2 m 2 m3 mim 2 mirri3 m\mz m 2 m3 / \m\m2 mimz m\mz , 



One may verify that the same expressions can be obtained by using (75). 



5 On controllability and observability of the physical systems con- 
sidered 



Theorem 5.1 For N = 2 the mechanical system £ is completely reachable 
and completely observable for any choice of mi, m 2 , C\ and hi. 

Proof. In connection with the square matrix F 2 of order 4 defined in (15), con- 
sider the characteristic polynomial ^(z) and the adjoint polynomial h 2 adj(a 2 ) g 2 , 
respectively defined by equalities (38) and (29) through (32). 

By Theorem 4.1 and Theorem 3.1, we respectively have 

4 , / i \ / 1 1 \ 2 7 t / \ / \ ZCi + fa 



P 2 (z) = z + (zci + fci)( 1 )z , h 2 adj(a 2 ) g 2 (z) 

mi m 2 mim 2 

Hence by replacing in P 2 (z) the unique root z\ = —ki/ci of h 2 adj(a 2 ) g 2 (z), we 
have P 2 {zi) = zf = (^j-) 4 7^ 0. Thus the characteristic polynomial and the ad- 
joint polynomial have no common root and Theorem 1.1 implies that £ is com- 
pletely reachable and completely observable; moreover, by assertions (E), (F) 
of Section 1.3, S is completely controllable and completely re-construictible 
too. 

Q.E.D. 

Theorem 5.2 Let N be any integer > 2. If 

ci/h = c 2 /k 2 = . . . = cjv-i/fcjv-i , (76) 



then 

(A) the mechanical system £ is completely reachable and completely observable 
for any choice of mi, m 2 , . . . m^. 

Proof. By Theorem 3.1 the roots of the adjoint polynomial (32) are Zi = 
—ki/ci, i — 1, . . . , N — 1. Let 

k = bi(z) = zci + h fa = -kj — + h , (77) 



17 



hence bi = if and only if z = z\. By hypothesis (76) we have 
/% = 0, i, j = l,2,...,N-l. 



(78) 



Hence the inductive definition (40)-(42) of Pn(z) yields 

P N ( Zj ) = zjP^zj) = z*P N _ 2 ( Zj ) = ... = (kj/ Cj ) 2N ? . (79) 

Thus any root Zj of the adjoint polynomial cannot be a root for the charac- 
teristic polynomial too. 

Now, Theorem 1.1 implies that £ is completely reachable and completely ob- 
servable; and by assertions (E), (F) of Section 1.3, £ is completely controllable 
and completely re-construictible too. 

Q.E.D. 

In particular, by the above theorem, for N = 3 the mechanical system E is 
completely reachable and completely observable for any choice of m 1? m 2 , 
7713, ci, ki, c 2 and k 2 provided that c\jk\ — c 2 /k 2 . 

The next two theorems respectively show that when the sequences {q} and 
{ki} are not proportional proposition (A) can be either true or not true. Hence, 
in particular, condition (76) is sufficient but not necessary for the validity of 
(A). 

Theorem 5.3 Proposition (A) is not true for some mechanical system E such 
that (76) does not hold. 

Proof. Let N = 3. We have 
i.e., 

Ps(z) = (V + ±) (V + b,— , 2 + z^\ + z^ U + h-) , (80) 
\ m 3 J \ m 2 m 1 J m 2 \ m 1 J 

Thus, if Z\ = is root of the adjoint polynomial, i.e. 

if &! = , then P 3 ( Zl ) = z\ + (3 12 Hz\ , with H=± + ± 

Hence ^3(21) = if and only if zj + f3 2 ±H = , i.e., 
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% + (k 2 - ^c 2 )H = 

C-y C\ 

This relation is linear in k 2 ; hence, if H, ki,c±, c 2 are given, then one determines 
the unique k 2 such that the last equality holds. With regard to these values of 
the material parameters it follows that Z\ is a root for both the characteristic 
polynomial and adjoint polynomial. 

Q.E.D. 

Theorem 5.4 Proposition (A) is true for some mechanical system £ such 
that (76) does not hold. 

Proof. Let iV = 3. By (80) we have 

Ps(*i) = U + ^) 4 + 4^4 = k 4\ k 4 + (k 2 ~ ^c 2 )H 23 ] , (81) 



m 3 l m 2 cf 



c\ Ci 



Ps(z 2 ) = 4 ( 4 + 4 + zt-^z\ = 4 \ k 4 + (*i - ^ Cl )H 12 , (82) 

V m 3 I m 2 c 2 - 



c 2 c 2 



where 



Ha = — + — (83) 



and 



Zi = -- (84) 



Ci 



are the roots of the adjoint polynomial. Hence 

P 3 ( Zl ) =0 & % + {k 2 - —c 2 )H 23 = (85) 
c i c i 



and 



Ps(z 2 ) = & % + (h - - Cl )H 12 = . (86) 

C 2 C 2 



Thus, given m 1; m 2 , m 3 , C\ and c 2 one can choose (hi, k 2 ) such that both 
-^3(^1) 7^ and P 3 (z 2 ) 7^ 0. Hence the adjoint polynomial and the characteristic 
polynomial have no common root. In connection with such a choice for the 
mechanical system E proposition (A) is true. 
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Q.E.D. 



6 On applications of the above results 

The system of equations Sjy studied here is a mathematical model suitable 
to describe various rectilinear physical systems as chains of three-dimensional 
bodies undergoing translational motions. 

By the results of Theorems 5.1-5.2 above, people involved in the design of 
such physical systems with N > 2 can choose the material constants in order 
to render the system completely controllable, observable, reachable and re- 
construictible; instead for N = 2 any such a system has such properties for 
any choice of the material constants. 

The key advantage of the theorems proved here for engineering applications is 
that, when control is a property useful for the system under design, one can 
design the system by choosing the material parameters in such a way as to 
have controllability. 

A further theoretical task will be to construct or to find the controlling input 
functions for the system. The present paper clarifies when such a problem has 
a solution. 

An example of a physical system, which can be modelled by T, N , is given 
by a chain of N pistons consecutively connected by springs and constrained 
inside a cylindrical cavity containing a fluid. Assuming that an external input 
force u{t) can act normally on the first piston, and taking as output the 
position of the latter piston, the system will be controllable and observable by 
the history of the latter provided the elastic springs be chosen proportional to 
the constants of viscous damping between consecutive pistons. Hence, (i) the 
system can reach any prescribed state by a suitable input force u(r) ; and (ii) 
every state of the system (that is, position and velocity of each piston) can be 
determined by observing the excursions of the latter piston. 

In particular, system S 2 can also be used to model the behaviour of car-wheel 
suspensions. In more detail, when a car travels along a bumpy road the wheel 
tyre copies roughness of the road surface. The wheel is thus driven up or down 
in the vertical direction along the z axis. In this case, the rectilinear system 
under consideration consists of four bodies: the spring, the shock absorber, 
the wheel and the quarter of the car body. The system excitation by the 
road surface and the gravitational attraction of the wheel and the quarter-car 
body form the system surroundings. The model is excited by the source of the 
velocity z Q resulting from the tyre copying the road bumps. 
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Figure 3 on the left shows a detail of the suspension of one of its wheels, and 
on the right shows a scheme of it. The point R z represents the z-position of 
the contact point between the car tyre and the road surface R when the car 
moves in the x-direction. The points W z and B z represent the z-motion of the 
wheel axis W and the body-suspension interaction point B, respectively. The 
constants k and c are the compliance and damping of the wheel tyre, mi is 
the mass of the complete wheel and m 2 is the quarter of the mass of the car 
body. The model is excited by the source of the velocity z Q resulting from 
the tyre copying the road bumps. Then by applying Newton's second law we 
obtain the differential equations of this physical model : 

d 2 

mi— z x = h(z 2 - Zl ) + Ci(i 2 - Zl ) + f{ xt , (87) 
d 2 

m 2 —z 2 = k 1 (z 1 -z 2 ) + c 1 (z 1 -z 2 ), (88) 



with 

/r* = k(z Q - Zl ) + c(z - ii) , (89) 

We see that the system considered is a particular S 2 system (see equations (7)- 
(9)). Hence by Theorem 5.1 it is observable and re-construictible for arbitrary 
choices of all the material parameters. 

If one wants to know how e.g. the wheel center W behaves in response to the 
external force road input, without a priori knowing the latter, he may observe 
the output past excursions y(s) := z 2 (s) of the body-suspension interaction 
point: they determine the actual state of the system, thus the actual position 
of W too. Thus he can go back to the external force road input which generates 
the vertical motion of the suspension. 
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k, k 2 
p OAA/WWXAA; ff OAA/WIA/IAA; Q 

Fig. 1. Physical model of system £2. 
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